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Abstract 

We study the combined effects of aging and links removal on epidemic dynamics in the Barabasi- 
Albert scale-free networks. The epidemic is described by a susceptible-infected-refractory (SIR) 
model. The aging effect of a node introduced at time ti is described by an aging factor of the form 
{t — ti)~^ in the probability of being connected to newly added nodes in a growing network under 
the preferential attachment scheme based on popularity of the existing nodes. SIR dynamics is 
studied in networks with a fraction 1 — p oi the links removed. Extensive numerical simulations 
reveal that there exists a threshold pc such that for p > pc, epidemic breaks out in the network. For 
p < Pc, only a local spread results. The dependence of pc on (3 is studied in detail. The function 
Pcil3) separates the space formed by /5 and p into regions corresponding to local and global spreads, 
respectively. 

PACS numbers: 89.75.-k, 89.75.Fb, 87.23.Ge, 87.19Xx 
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I. INTRODUCTION 



It has been discovered in recent years that many real network systems, while showing 
different levels of complexity of their own, possess novel common structural or topological 
properties. These networks include, for example, the Internet World- Wide- Web 2|, 
scientific citations i, cells Q, the web of collaborations among actors and actresses j5|, 
and the web of human sexual contacts [6J . Typically, networks are characterized by nodes 
and links. The nodes may represent websites on the world- wide- web, a paper in scientific 
journal, or an actor, etc., depending on the network concerned. The links represent the 
interaction between the nodes. The interaction may represent the links from a webpage to 
another, references cited in a journal to papers in the literature, collaboration among actors 
and actresses, etc., again depending on the system concerned. These structures are found 
to show the small- world effect, i.e., a node can reach any randomly chosen node through 
only a few links; the high-clustering effect, i.e., nodes connected to a chosen node have high 
probability of having links among themselves; and a well-defined degree distribution |2,y|- 
These universal features observed in real systems have led to the intensive studies of complex 
networks in recent years. Physicists have made major contributions to the understanding of 
the physics of networks. On the structural properties, physicists have proposed a number 
of models in which the observed statistical features in real networks are reproduced. The 
most representative of these models are the small-world model proposed by Watts and 
Strogatz ; and the scale- free growing network proposed by Barabasi and Albert (BA) jl^ . 
The latter model gives a better description of the scale-free or broad-scale distributions of 
degrees found in many real-life networks and hence has been extensively studied recently. 
While structural properties are important in the geometrical description of networks, many 
applications of the physics of networks rely on our understanding of dynamical processes on 
networks. Dynamical processes, such as percolation, searchability problems, and especially 
the spread of diseases on complex networks have been investigated only recently. It has been 
fouud that the BA networks are highly susceptible to large scale epidemics Q. 

In the present work, we study how aging phenomena and links removal combined 
may affect the spreading of a disease or rumor in BA growing networks. Aging refers to 
situations in which the older nodes in a network become increasingly unattractive to the 
newly added nodes. We use a susceptible-infected-refractory (SIR) model for the dynamics 
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of epidemics. Using extensive numerical simulations, we study the SIR dynamics on BA 
networks with a fraction of 1 — p of the links removed or turned ineffective. It is found that 
a threshold pdP) exists, where f3 is an aging parameter. For p > Pc, 9. global spread results; 
while for p < pc, only a local spread results. For the full BA network, i.e., p = 1, global 
outbreak can take place only if the value of f3 is smaller than a critical value. 



II. MODEL 



The underlying network can be set up as follows. Initially there are tuq fully connected 
nodes. We take mo = 5 in the present work. At each time step, a new node is introduced. 
Each newly added node establishes m (< mo) outgoing links to existing nodes. Each new 
link introduced at time step t has a probability 
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to be connected to an existing node i with degree fc,, i.e., having ki links. Here is the 
time at which node i was added to the network and (3 is an aging parameter characterizing 
how rapidly an existing node becomes unattractive to a newly added node. Eq.((H) states 
that newly established links are preferentially attached to the younger and highly connected 
existing nodes in the system. For /5 = 0, there is no aging effect and the model reduces 
to the BA model 2. The aging factor thus models the situation where older persons in 
a society tend to isolate themselves or being isolated by the younger ones and lose their 
influence in the younger persons. 

As an undiluted BA network is highly susceptible to large scale global spread of diseases, 
it is interesting to study the effects of removing a fraction of links in a fully grown BA 
network. These diluted networks can be constructed as follows. Having set up an aging 
network of nodes, each link has a probability 1 — p of being removed. Hence, on the 
average, a fraction p of the links in a fully grown aging network consisting of N nodes are 
kept. The removed links in a diluted network can be regarded as ineffective links through 
which a disease or a rumor cannot spread. In the present work, we study aging networks of 
size up to = 50, 000 nodes. 

To study epidemiological processes on aging and diluted networks, we use the three- 
state Susceptible- Infected- Refractory (SIR) model [l^. The model is a standard model 
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for studying epidemics and the spread of a rumor in a connected population. Initially, all 
nodes are in the S(susceptible)-state, and one node is randomly chosen to be infected, i.e., 
in the I-state. The SIR dynamics proceeds as follows. At a time step t, a node, say i, is 
randomly chosen among all the infected (1) nodes. A neighbor or friend j is then selected 
randomly among all the neighbors of node i, i.e., those with a link connected to i. If node 
j is susceptible, it becomes infected and the chosen node i remains in the I-state; otherwise 
(i.e., node j is either I or R) the state of node j remains unchanged and the chosen node 
i becomes refractory (R) at the end of the time step. Within the context of a spread of a 
rumor, the I-state nodes refer to persons who want to spread the rumor, the S-state nodes 
arc persons who have not heard the rumor, and the R-state nodes refer to persons who lost 
interest in spreading the rumor after knowing it. As time evolves, the number of R- nodes 
(S-nodes) increases (decreases); while the number of I- nodes increases initially and then 
eventually drops to zero. The number of R-nodes at the end of the dynamical evolution is 
denoted by Nr. We are interested in how Nr/N varies as the fraction of effective links p 
and the aging parameter /3. 



III. RESULTS AND DISCUSSION 



We performed extensive numerical simulations to explore the combined effects of aging 
and links removal on the SIR dynamics. In all the simulations, each data point is obtained 
by averaging over 100 realizations of the network structure and 100 different initially infected 
nodes for each realization. In Fig.l we show the mean fraction of R-nodes, r =< Nr > /N, 
on networks of different sizes up to = 50, 000 on a log-log plot for different values of the 
aging parameter (a) /3 = 0, (b) /3 = 0.5, (c) j3 = 1.0, and (d) f3 = 1.5. For each value 
of (3, we carried out simulations for different levels of dilution p. For global spread, we 
expect that the number of refractory sites Nr to be directly proportional to the network 
size A^ and hence the fraction of refractory nodes r becomes independent of N. For local 
spreads, Nr is finite for a sufficiently large network and hence r decreases as 1/A^. We 
thus search for a threshold Pc that separate these two behavior of r for a given value of /3. 
The results in Fig. 1 show that r = Nr/N docs exhibit different behavior as p decreases. 
Precise determination of Pc is quite difficult, due to finite size effect and the inhomogeneous 
nature of the problem resulting from the random removal of links from a originally randomly 
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established BA network. To proceed, we study the slope of the lines in Fig.l as p decreases. 
Typical results for (3 = 0, 0.5, and 1.0 are shown in Fig.2. Notice that, for given value of 
/3, the slope starts to drop from zero, i.e., Nji ~ for global spread, for some value of p. 
Notice that there is a transitional range of p for which N^i neither scales linearly with 
nor independent of A^. We adopt the value of p at which the slope drops to the value of 
—0.1 as Pc{P)- According to this criteria, we found pc{P = 0) = 0.2, pc{P = 0.5) = 0.25, 
PciP = 1) = 0.3, and pdP = 1.5) = 0.35 for the values of /5 in Fig.l and Fig.2. 

Carrying out the calculations for different values of the aging parameter /3, we obtained 
the dependence of the threshold pdP) on the aging parameter p. The results of pdP) are 
shown in Fig. 3 as the phase boundary in a phase space formed by the aging parameter (3 
and fraction of effective links p. For small jS {0 < (3 < 1.25), Pc increases gradually as (3 
increases. For (3 > 1.25, pc increases more sensitively with (3. As (3 approaches 2, pc rapidly 
increases and takes on a value close to 1, indicating that a rapidly aging network has the 
slightly older nodes effectively isolated and hence only local spread of a disease or rumor 
is possible, regardless of whether links are further removed. The phase boundary Pc{(3) 
separate the phase space into two regions, one corresponds to local spread of a disease and 
another corresponds to a global spread of a disease. Note that for (3 > Pc, where Pc ~ 2.0, 
Pc = ^ implying that aging effect alone is sufficient to keep a disease or rumor from spreading 
globally. 

The effect of aging can be understood qualitatively as follows. A BA network without 
aging, i.e. /? = 0, is highly resilient to random damages in the form of link removal Q]. The 
results in Fig. 1(a) show that one needs to remove more than 80% of the links to prevent 
global spread of an epidemic. This feature is similar to previous studies of site percolation 
15[, where it was shown that for networks with a broad-scale degree distribution, one 
needs to remove a large fraction of nodes (about 99%) before the network falls apart. The 
sensitivity of pc to the aging parameter P is related to the form of the degree distribution 
as P increases. For small values of P {P < 1), the degree distribution still carries scale-free 
features that lead to a small pc and the weak dependence of pc on p. For larger values of P, 
the older nodes lose their attractiveness rapidly giving rise to a narrow degree distribution. 
The scale-free nature of the network is lost and the degree distribution becomes increasingly 
dominated by an exponential form rather than a power law. In this case, local spread is 
more probable unless the network carries a larger fraction of effective links, giving rise to 
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the sensitive dependence of pc on (3 0, Ol- For rapidly aging networks, the small world 
effect vanishes and the "regular lattice" behavior dominates. The hig h clustering feature of 
a rapidly aging network results in only local spread of a diseases jl6 |. 

In summary, we have studied the combined effects of aging and links removal on SIR 
dynamics in the Barabasi-Albert scale-free networks. A threshold Pc exists for the fraction 
of effective links in a network below which only local spreads of disease or rumor take place. 
For small aging effect, pc is low. As aging effect increases, pc increases. For rapidly aging 
network, Pc=l implying that only local spread results, regardless of the number of effective 
links in the network. 
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FIGURE CAPTIONS 



Fig. 1: The fraction of refractory nodes r = Nr/N as a function of the size of network N for 
different fractions of effective hnks in BA networks on a log-log scale. The aging parameter 
is (a) /3 = 0, (b) /3 = 0.5, (c) f3 = 1.0 and (d) /3 = 1.5. Note that for global spreads, Nr ^ N 
and a horizontal line results. 

Fig. 2: The slope of logr against log A as a function of fraction p of effective hnks in a 
network. Note the change from vanishing slope to negative slopes as p decreases, signifying 
a change from global to local spread of an epidemic. 

Fig. 3: The thresholds pdP) constitute a boundary in the phase space formed by the aging 

parameter f3 and fraction of effective links p. In the region below (above) the boundary, 
local (global) spreads occur. 
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